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This  paper  studies  the  transient  response  of  the  output  voltages  of  a  Ballard-Mark-V  35-cell  5  kW  pro¬ 
ton  exchange  membrane  fuel  cell  (PEMFC)  stack  with  power  conversion  for  applications  in  autonomous 
underwater  vehicles  (AUVs)  under  load  changes.  Four  types  of  pulse-width  modulated  (PWM)  dc-dc 
power  converters  are  employed  to  connect  to  the  studied  fuel  cell  in  series  for  converting  the  unregu¬ 
lated  fuel  cell  stack  voltage  into  the  desired  voltage  levels.  The  fuel  cell  model  in  this  paper  consists  of  the 
double-layer  charging  effect,  gases  diffusion  in  the  electrodes,  and  the  thermodynamic  characteristic; 
PWM  dc-dc  converters  are  assumed  to  operate  in  continuous-conduction  mode  with  a  voltage-mode 
control  compensator.  The  models  of  the  study’s  fuel  cell  and  PWM  dc-dc  converters  have  been  imple¬ 
mented  in  a  Matlab/SIMULINKTM  environment.  The  results  show  that  the  output  voltages  of  the  studied 
PEMFC  connected  with  PWM  dc-dc  converters  during  a  load  change  are  stable.  Moreover,  the  model 
can  predict  the  transient  response  of  hydrogen/oxygen  out  flow  rates  and  cathode  and  anode  channel 
temperatures/pressures  under  sudden  change  in  load  current. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Over  the  past  few  years,  hydrogen  production  and  fuel  cell  tech¬ 
nologies  have  grown  dramatically.  Proton  exchange  membrane  fuel 
cells  (PEMFC)  are  one  of  the  most  promising  fuel  cell  technologies 
that  can  be  a  candidate  for  both  stationary  and  automotive  appli¬ 
cations  [1].  The  traditional  energy  source  of  underwater  vehicles 
is  based  on  the  heat  engine  process  of  the  Carnot  cycle.  It  causes 
not  only  low  conversion  efficiency  but  also  pollutions.  Recently, 
an  unmanned  untethered  submersible  AUV-HM1  testbed  devel¬ 
oped  by  National  Taiwan  University  has  been  powered  by  lead-acid 
batteries.  However,  the  limited  operating  range  detracts  from  the 
usages  of  them.  In  order  to  design  a  long-range  underwater  vehi¬ 
cle,  the  larger  capacity  of  the  lead-acid  battery  is  needed  but  this 
will  increase  the  weight  and  occupy  most  of  the  volume  of  the 
underwater  vehicle.  Thus,  fuel  cells  coupled  with  an  alternative 
reactant  storage  system  may  have  tried  to  overcome  these  limita¬ 
tions  while  maintaining  other  benefits  of  an  electrochemical  power 
source  [2-6]. 

Generally,  the  PEMFC  has  multiple  advantages.  It  offers  a  quiet 
low  temperature  operating  condition  with  only  producing  a  clean 
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byproduct— pure  water  and  is  stackable  to  provide  a  range  of  space 
options.  Its  disadvantage  is  that  pure  hydrogen  and  oxygen  reac¬ 
tants  must  be  carried.  Thus,  fuel  storage  is  a  key  component  for 
applying  the  PEMFC  to  underwater  vehicles.  For  the  PEMFC,  pure 
hydrogen  must  be  stored.  Hydrogen  can  be  stored  as  a  gas,  liquid, 
or  in  a  hydride.  Out  of  the  three  options  researched,  the  hydride 
method  is  chosen  for  its  proven  safety  and  reliability  advantages 
[4]. 

The  purpose  of  this  paper  is  to  study  the  performance  of  the  com¬ 
bined  PEMFC  and  pulse-width  modulated  (PWM)  dc-dc  converters 
subject  to  a  load  change  in  autonomous  underwater  vehicles 
(AUVs).  Although  [7]  has  already  described  the  PEMFC  performance 
in  the  presence  of  transient  effects,  it  did  not  use  power  convert¬ 
ers  to  regulate  the  output  voltage  of  the  fuel-cell  stack.  In  addition, 
[8]  has  built  a  nonlinear  model  of  buck  converters,  and  in  [9-11] 
have  mentioned  to  the  use  of  boost  converters  to  regulate  the  out¬ 
put  voltage  of  the  fuel-cell  stack.  Nevertheless,  they  have  not  yet 
shown  the  phenomenon  that  the  load  variations  may  affect  the 
output  voltage  of  the  power  converter. 

2.  Mathematical  model  of  the  PEMFC 

This  section  describes  a  mathematical  model  for  investigating 
the  mass  and  energy  dynamic  interactions  as  well  as  the  electro¬ 
chemical  reaction  within  a  PEMFC  [12,13].  As  seen  from  Fig.  1,  the 
hydrogen  enters  the  anode  channel,  and  diffuses  through  the  anode 
diffusion  layer  toward  the  active  catalyst  layer,  where  the  hydrogen 
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Fig.  1.  Schematic  diagram  of  a  PEMFC. 


dissociates  into  protons  and  electrons.  The  protons  then  travel  from 
the  anode  catalyst  to  the  cathode  catalyst  through  a  membrane.  In 
the  meantime,  oxygen  enters  the  cathode  channel,  part  of  which 
diffuse  through  the  cathode  diffusion  layer  to  the  active  catalyst, 
where  the  oxygen  dissociates  and  combines  with  protons  and  elec¬ 
trons  to  form  water,  while  heat  and  current  are  generated.  The  heat 
produced  during  the  reduction  reaction  process  will  warm  up  the 
fuel  cell  body  that  will  affect  the  mass  diffusion  within  the  anode 
diffusion  layer,  cathode  diffusion  layer  and  membrane.  Typically, 
the  heat  transfers  to  ambient  and  gas  flow  channels  through  con¬ 
vection  and  radiation.  The  resulting  gas  temperature  variation  will 
change  the  pressure  and  flow  rate  within  the  channels  to  affect  the 
gas  diffusion  within  the  diffusion  layers.  All  these  dynamic  inter¬ 
actions  will  have  a  direct  effect  on  the  fuel  cell  output  voltage  and 
current. 


cell.  In  addition,  the  physical  properties  associated  with  the  control 
volumes  are  assumed  to  be  uniform. 

2.1  A.  Anode  and  cathode  channels 

Generally,  air  or  oxygen  is  fed  to  the  cathodic  compartment 

while  hydrogen  is  fed  to  the  anodic  one  and  the  electrolyte  per¬ 
forms  both  the  functions  of  transferring  H+  from  the  anode  to  the 
cathode  and  reactant  separation.  Note  that  pure  hydrogen  and  oxy¬ 
gen  are  served  as  the  fuel  and  oxidant  gases  for  underwater  appli¬ 
cations,  respectively.  Thus,  the  first  term  in  (1)  becomes  the  mass 
change  rate  of  the  respective  gases  for  the  control  volumes  corre¬ 
sponding  to  the  anode  and  cathode  channels.  The  second  term  in 
(1 )  is  the  mass  flow  rate  through  the  surface  of  the  control  volumes. 

For  the  anode  channel,  the  hydrogen  mass  change  rate  is  written 
to  be  [12,13]: 


2.1.  Relationship  between  mass  and  energy  transfer  using  control 
volume  approach 


c.v. 


diriH2 

dt 


(3) 


For  an  arbitrary  control  volume,  a  continuity  equation  in  an 
integral  form  is  given  as  [14]: 


and  the  mass  flow  rate  through  the  surface  of  the  control  volumes 
is  expressed  by 


d_ 

dt 


p&V+  §p{V 

C.S. 


h)dA  =  0 


c.v. 


(1) 


and  an  energy  equation  in  integral  form  is  expressed  by 

C.V. 


where  p  is  the  material  density  within  the  control  volume,  V  is  the 
velocity  vector  of  flow  fluid,  n  is  the  unit  vector  normal  to  the  con¬ 
trol  surface  toward  the  outside,  et  is  the  specific  internal  energy, 
h  is  the  specific  enthalpy,  Q.  is  the  heat  added  to  the  control  vol¬ 
ume,  and  W  is  the  work  done  on  the  control  volume.  Since  the  gas 
flow  velocity  is  usually  low  under  fuel  cell  operating  conditions,  the 
kinetic  and  potential  energies  of  the  gases  can  be  neglected  in  (2). 
Note  that  the  anode  channel  control  volume  is  defined  to  be  the 
lumped  volume  from  the  anode  inlet  to  the  anode  outlet.  Similarly, 
the  cathode  channel  control  volume  is  defined  at  the  cathode  side. 
The  membrane  electrode  assembly  (MEA)  and  all  other  metal  struc¬ 
tures  involved  in  the  fuel  cell  will  form  the  fuel  cell  body  control 
volume.  As  for  considering  a  fuel  cell  stack,  the  control  volumes 
of  the  anode  channel,  cathode  channel  and  fuel  cell  body  will  be 
the  summation  of  corresponding  control  volumes  of  each  single 


-fsPH2{V  n)dA  =  WHl-m  -  Wh2, out  -  WH2,used  (4) 

where  in  and  WH2,0ut  are  the  inlet  and  outlet  mass  flow  rates, 
respectively.  They  are  as  functions  of  the  hydrogen  source  pressure, 
the  pressure  within  the  anode  control  volume  and  ambient  pres¬ 
sure.  With  the  macroscopic  level  dynamic  behavior  of  the  fuel  cell, 
the  inlet  and  outlet  mass  flow  rates  can  be  simplified  by  using  the 
traditional  nozzle  flow  rate  equation  which  are  written  to  be  [12]: 


^H2,in  —  ^an,in(Pan,s  “Pan) 

(5) 

k^H2,out  =  ^an,out(Pan  “Patm) 

(6) 

where  /<an  in  and  /<an,out  are  the  inlet  and  outlet  mass  flow  rate  coef¬ 
ficients;  pan,s,  Pan  and  patm  are  the  hydrogen  source  pressure,  the 
pressure  within  the  anode  control  volume  and  ambient  pressure, 
respectively.  In  addition,  the  consumed  mass  flow  rate  of  hydrogen 
is 

Moused  =MH2%r  (?) 

where  /  denotes  the  current,  Ncell  is  the  number  of  cells  in  a 
stack;  Fand  MH2  are  Faraday’s  constant  and  hydrogen’s  mole  mass, 
respectively.  Thus,  the  hydrogen  mass  transport  rate  through  the 
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anode  channel  control  surface  becomes: 

—  ^an,in(Pan,s  —  Pan)  —  ^an,out(Pan  —  Patm)  —  ^U2  2f  (^) 

Furthermore,  the  energy  equation  of  (2)  can  be  applied  to  the  anode 
in  a  similar  manner.  The  first  term  on  the  left-hand  side  of  (2)  repre¬ 
sents  the  internal  energy  change  rate  of  hydrogen  within  the  anode 
channel  control  volume  which  is  given  as  [12,13]: 


_9 

dt 


Pn7et  dV 


c.v. 


d(mH2Cy;H2^an) 

dt 


(9) 


where  cv^2  is  the  specific  heat  of  hydrogen  at  constant  volume  and 
Tan  is  the  hydrogen  temperature  within  the  anode  channel  control 
volume.  The  second  term  on  the  left-hand  side  of  (2)  is  the  heat 
transfer  rate  due  to  the  mass  transport  through  the  control  surface 
which  is  given  as: 

— Ph2 ^H2 iy  ’  n)dA  =  V^H2,in^H2,in 


^H2,out^H2,out  ^H2,  used  ^H2,  used 


(10) 


d(mo2  Cv,02  Tea  ) 

dt 


kconv,  ca  (  T)ody  Tc  a)  +  ^ca,in(Pca,s  Pca)Cp,02  (Tea  ,in  To) 

~^ca,out(Pca  —  Patm)Cp,02  (TCa  —  To)  —  ^02Cp, 02(TCa  —  To) 


(19) 


2.2.2.  Cell  body 

The  continuity  equation  for  the  cell  body  control  volume  is  only 
concerned  with  the  species’  conservation  of  electrochemical  reac¬ 
tion.  That  is,  water  is  formed  through  the  electrochemical  reaction 
by  consuming  hydrogen  and  oxygen.  Since  most  parts  of  the  cell 
body  are  composed  of  metal  structures,  gas  diffusion  and  water 
produced  have  little  effect  on  the  total  mass  change  of  cell  body. 
Thus,  we  will  have  [12,13]: 


_9 

dt 


Pbody  dV  —  0 


C.V. 


(20) 


and  the  mass  transport  rate  through  the  cell  body  control  surface 
(Membrane  Electrode  Assembly,  MEA)  is 

-fsPbody(V-n)dA  =  WH 

2,used  T  Wq2  ?  used  VVh2  o  (21) 


where  ftH2,in  and  fiH2,out  are  the  inlet  and  outlet  specific  enthalpy 
hydrogen,  respectively;  hH2 ?used  is  the  specific  enthalpy  of  con¬ 
sumed  hydrogen.  Since  the  temperature  of  a  PEMFC  body  is 
relatively  moderate,  the  radiation  effect  can  be  neglected.  Based 
on  Newton’s  cooling  law  for  governing  the  convection  heat  trans¬ 
fer  mechanism,  the  convection  heat  transfer  from  the  cell  body  to 
the  anode  channel  can  be  described  as  [12,13]: 

—  Qconv,an  —  hconv,anAconv,an  AT  =  kconv,  anybody  —  Tan)  (11) 

where  /iCOnv>anAconv,an  =  kconv,an»hConv,an  is  the  convection  heat  trans¬ 
fer  coefficient  at  the  anode;  Aconv,an  is  the  surface  area  of  hest 
convection  at  the  anode.  Note  that  the  work  done  on  the  control 
volume  (i.e.,  dW/dt)  disappears  since  the  viscous  normal  stresses 
at  the  inlet,  outlet  and  electrode  interface  of  mass  diffusion  are 
neglected.  Thus,  the  energy  equation  for  the  anode  becomes: 


d(mH2  c„,H2Tan) 

dt 


fcconv,an(Tbody  Tan)  +  /Can,in(Pan,s  Pan)Cp,H2  (Tan, in  To) 

— ^an,out(Pan  —  Patm)Cp,H2  (Tan  —  To)  —  ^  ^H2  Cp, H2  ( Tan  —  To) 

(12) 

where  Tan  in  and  T0  are  the  flow-in  hydrogen  temperature  and  the 
reference  temperature  (i.e.,  298.1 5  K),  respectively.  cp?H  is  the  aver¬ 
age  specific  heat  of  hydrogen. 

For  the  cathode  channel,  similar  equations  can  be  derived  by 
replacing  the  relevant  subscripts  by  “ca”  (cathode)  and  “02”  (oxy¬ 
gen).  The  following  equations  are  parallel  to  those  given  in  (3)  and 
(4)  and  (8)— (12): 


8 


dt  J  1  1 
c.v 


Po,  dV  = 


dm 


o2 


dt 


fcPo2(V  '  ft)dA  —  Wo2,in  Wo2,out  Wo2,used 


c.s 

dm 


(13) 


(14) 


—  ^ca,in(Pca,s  —  Pea)  —  kca,out(Pca  —  Patm)  —  Mo  ^  (15) 


dt 

dt  J  P°2 

C.V. 


4F 


Po2et-dV  = 


d(m02Cp?o2^ca) 

dt 


(16) 


-§Po2ho2{y  •  n)  •  dA  =  Wo2,in^o2,i 
-Wo2 

,out^02,out  i^02,  used  ^02,  used 

da 


dt 


—  Qconv,ca  —  kConv,ca(l"body  —  fca) 


(17) 

(18) 


Note  that  the  above  equation  includes  the  gas  diffusion 
through  the  electrodes  and  water  formed  by  electrochemical 
reactions. 

Now  let  us  consider  the  energy  equation  for  the  cell  body  control 
volume.  The  internal  energy  change  rate  of  cell  body  is  given  as 
[12,13]: 


d_ 

dt 


dr, 


Ilf  uu 

Pbody^t  Ov  =  mbodyCp  body 


body 


(22) 


c.v. 


and  the  heat  transfer  rate  due  to  the  mass  transport  through  the 
control  surface  is  given  as: 


c.s. 


Pbody^body(^  '  tl)dA  —  Wh2, used ^H2, used 


“T  ^02, used ^02, used  ^H20^H20 


(23) 


where  cpbody  is  the  average  specific  heat  of  the  cell  body.  In 
addition,  several  heat  transfer  effects  are  involved  in  the  term 
dQ/dt,  which  include  the  heat  added  to  the  cell  body  through 
the  electrochemical  reaction  as  well  as  the  heat  transferred  to 
the  ambient  and  the  channels  via  convection.  Thus,  we  have 
[12,13]: 


dQ 

—  Qconv,an  +  Qconv,ca  +  Qconv,amb 
+  A  Hrj  #  PH2(V-n)cL4 

electrode 


(24) 


where  Qconv.an  is  the  heat  convection  between  the  cell  body  and  the 
anode  channel,  Qconv.ca  is  the  heat  convection  between  the  cell  body 
and  the  cathode  channel,  Qc0nv,amb  is  doe  heat  convection  between 
the  cell  body  and  the  ambient,  and  A HR  T  is  the  lower  heating  value 
of  the  fuel.  All  the  heat  convection  terms  here  are  assumed  to  com¬ 
ply  with  Newton’s  cooling  law.  For  the  cell  body,  the  term  dW/dt  is 
the  electric  power  converted  from  electrochemical  reaction,  i.e.: 


=  -V1  I 
dt  stack 


(25) 


where  V'  k  denotes  the  stack  voltage.  Thus,  the  heat  added  to  and 
the  work  done  on  the  cell  body  control  volume  is  then  characterized 
as: 


dQ.  dW 

—  Kconv,anUan  —  ^body)  +  onv,ca(4ca  —  JbodyJ 

+  kconv,amb(Tamb  "  Tbody)  +  AHRj  ^Mh2  -  V'^l  (26) 
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Note  that  the  first  three  terms  at  right-hand  side  are  the  effects 
of  the  convection  heat  transfer  from  the  surrounding  structures, 
the  fourth  term  on  the  right-hand  side  is  the  heat  produced  by  the 
electrochemical,  and  the  last  term  is  the  converted  electric  power. 
Finally,  the  energy  equation  for  the  cell  body  will  become 


As  the  partial  pressure  of  water  vapor  in  the  cathode  is  fixed  at  the 
saturation  level,  the  partial  pressure  of  02  at  the  interface  is 

Pho 

Po2  =(1-Xh2o)^  (34) 

Xh2o 


^body^p,body 


d^body 

dt 


A Hrj^-MH2  - 
+WH,  ,used^p,H2  (Tan 


^staclJ  +  kconv,an(7an  —  Tbody)  +  kconv,ca(Tca  —  ?body)  +  ^conv,amb(^amb  —  ^body)  (27) 
—  ^body)  T  ^02,used^p,02  (^ca  —  ^body)  —  ^H2oCp,H2o(7ca  ~~  ^body) 


2.2.  Reactant  gases  diffusion  in  the  electrodes 

To  simply  the  mathematical  model  of  the  PEMFC,  the  following 
assumptions  are  made  in  this  paper  [15-21]: 


Also,  the  effective  cathode  water  vapor  pressure  equals  the  satu¬ 
rated  vapor  pressure  based  on  assumption  4.  Thus,  the  effective 
partial  pressure  of  02  by  combining  (33)  and  (34)  yields: 


(1)  One-dimensional  treatment. 

(2)  Ideal  and  uniformly  distributed  gases. 

(3)  Constant  pressures  in  the  fuel-cell  gas  flow  channels. 

(4)  The  fuel  is  humidified  H2  and  the  oxidant  is  humidified  02. 
Assume  the  effective  anode  water  vapor  pressure  is  50%  of 
the  saturated  vapor  pressure  while  the  effective  cathode  water 
pressure  is  100%. 

(5)  The  reaction  gases  are  pure  hydrogen  and  oxygen. 

(6)  Parameters  for  individual  cells  can  be  lumped  together  to  rep¬ 
resent  a  fuel-cell  stack. 

(7)  As  mentioned  in  [20],  the  saturated  vapor  pressure  pjjat0  is  a 
function  of  the  operating  temperature  of  the  PEMFC  and  is  writ¬ 
ten  to  be 


log(p^0)  =  2.95  X  10-2(Tbody  -  273.15)  -  9.18  x  10“5 

(Tbody  -  273. 15)2  +  1.44  x  l(T7(rbody  -  273. 15)3  -  2.18  (28) 

In  order  to  calculate  the  PEMFC  output  voltage,  the  effective  partial 
pressures  of  H2  and  02  need  to  be  determined.  In  a  gas  mixture 
consisting  of  N  species,  the  diffusion  of  component  i  through  the 
porous  electrodes  can  be  described  using  Stefan-Maxwell  equa¬ 
tions  as  written  to  be  [21  ]: 


VX/  = 


XM),  i  =  1,2, 


n 


(29) 


In  the  one-dimensional  transport  process  along  the  X-axis  as  shown 
in  Fig.  1,  the  effective  mole  fraction  of  water  at  the  anode  catalyst 
interface  is  given  as: 


X, 


h2o 


=  X“‘  exp 


^^body/^ac 

2FPan^H20,H2 


(30) 


Since  Xh2  +XH2o  =  1,  the  effective  partial  pressure  of  H2  at  the 
interface  is 

PSHO 

PH2=^-XH2o)tPX  (31) 

^h2o 

According  to  assumption  4  mentioned  above,  the  effective  anode 
water  vapor  pressure  is  50%  of  the  saturated  vapor  pressure.  Thus, 
the  effective  partial  pressure  of  H2  by  combining  (30)  and  (31) 
yields: 


Ph2  —  0.5^ 


x£aonne'  exp((RTbodyi4c)/2FPanDH2o,H2 ) 


-  1 


(32) 


Similar  to  the  analysis  for  the  anode,  the  effective  mole  fraction  of 
water  at  the  cathode  catalyst  interface  can  be  found  as: 


V  ^channel 

Xh2o=xho  exp 


ftTbod yjt 


cc 


4FPcaDH2o,02 


(33) 


Po2  -  P h2o 


1 

Xcharmd  exp((RrbodyJ€Cc  )/4 FPcaD  h20,02) 


1 


(35) 


Using  the  appropriate  critical  values  and  constants  for  a 
water/oxygen  mixture,  (32)  and  (35)  finally  becomes: 


pw2  —  0.5p^o 


X“exp((1.653J)/rio3d3y4) 


(36) 


Po2  -  Ph9o 


1 

XchaonneleXp((4.192;)/T1]o3d3y4) 


1 


(37) 


where  (36)  and  (37)  will  be  used  in  the  Nernst  equation  to  find  the 
fuel-cell  output  voltage. 


2.3.  Output  voltage  of  the  fuel-cell  stack 

The  total  output  voltage  of  a  fuel-cell  stack  is  a  function  of 
stack  current,  cathode  pressure,  reactant  partial  pressures,  fuel  cell 
temperature,  membrane  humidity,  etc.  Since  the  fuel-cell  stack  is 
comprised  of  multiple  fuel  cells  connected  in  series,  the  stack  volt¬ 
age  Vstack  is  calculated  by  multiplying  the  cell  voltage  Vcelj  by  a 
number  of  cells  Nce\\  in  the  stack,  which  is  written  as: 

^stack  =  ^cell  *  ^cell  (38) 

The  output  voltage  of  a  single  fuel  cell  is  defined  by  the  sum  of  three 
voltage  terms:  the  thermodynamic  potential  Ene mst,  the  activation 
overvoltage  Vact,  and  the  ohmic  overvoltage  Vohm  as  repeated  below 
[15-19]: 

Vcell  =  ^nernst  —  i^act  —  i^ohm  (39) 

where 

^nernst  =  Eq  —  0.85  x  10  ( T^body  —  298.15) 

+  RfpLMpH2-(p02)0-5}  (40) 


^act  =  -{§1  +  §2^body  +  §3^body  Kn(Q}2  )1  +  ?4?body^n(0} 

§2  =  0.00286 +  0.0002£n(Acell)  + 4.3  x  10“5€n(CH2) 


v  /  498  \ 

C02  =  1.97  xl0-7po2  exp  - - 

VTbody/ 


—77 

Ch2  =  9.174  x  10~7Ph2  exP 


body 


^ohm  —  J  •  ^int 

r  mpm  tr 


^int  — 


mem  Lmem 


aceii 


(41) 

(42) 

(43) 

(44) 

(45) 

(46) 
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181. 6[1  +  0.03J  +  0.062(rbody/303)2/2  5] 
rmem  [A-0.634-3J]exp[4.18((rbody-303)/Tbody)]  1 

Note  that  the  two  overvoltage  terms  in  (39)  are  both  negative 
and  represent  nonideal  fuel-cell  reaction  conditions.  From  an 
electrical/power-electronics  standpoint,  a  fuel  cell  is  almost  a  non¬ 
linear  current-dependent  voltage  source.  While  this  functional 
relationship  varies  somewhat  with  different  fuel  cells  (e.g.,  solid 
oxide  fuel  cell,  SOFC),  [22-24]  have  outlined  more  detailed  work 
on  planar  SOFC. 

2.4.  Double-layer  charging  effect 


The  double  layer  charging  capacitance  may  affect  the  activation 
overvoltage  term  in  (39).  This  results  in  a  first-order  model  where  a 
sudden  change  in  cell  current  is  associated  with  a  slower  variation 
in  output  voltage  [17,19].  An  electrical  circuit  model  is  generally 
used  to  illustrate  this  capacitive  effect  [17].  For  a  given  cell  current, 
the  activation  voltage  drop  can  be  considered  across  an  activation 
resistance  Rac t.  Since  the  activation  voltage  term  l/act  introduced  in 
(41 )  is  negative,  the  derivation  of  a  positive  term  Vd  is  required.  The 
first-order  dynamics  of  the  activation  overvoltage  are  described  by 
the  equations  as  follows  [17,19]: 

Ract  =  ^  (48) 

=  _L  _  _^L_  (49) 

dt  Cdl  RactQn  1  ; 

If  the  double-layer  charging  effect  is  integrated  into  the  modeling 
by  using  Vd  instead  of  Vact,  to  calculate  V'  k,  the  fuel-cell  output 
voltage  now  turns  out  to  be: 

=  Ncell(£nernst  -  Vd  -  Vohm )  (50) 

Note  that  the  dynamical  model  for  PEMFC  electrochemistry 
described  in  this  section  is  based  on  a  generalized  empirical  set  of 
equations.  Thus,  by  identifying  appropriate  parameters  (i.e.,  num¬ 
ber  of  cells,  cell  area,  membrane  thickness,  etc.),  fuel  cells  with 
varying  sizes,  structures,  and  materials  can  be  modeled,  and  the 
relevant  dynamics  can  be  studied  [17,19].  Table  1  lists  the  related 
parameters  used  to  model  the  PEMFC  in  this  paper. 


3.  Design  of  DC/DC  power  converters 

For  the  designed  AUV,  there  are  four  desired  voltage  levels 
including  5-,  12-,  24-  and  48-\/dc  in  response  to  changes  in  the 
output  load  and  the  PEMFC  output  voltage.  To  achieve  this,  a  nega¬ 
tive  feedback  control  system  is  used  where  the  converter  output  is 
compared  with  its  reference  value.  The  compensator  produces  the 
control  voltage,  which  is  used  to  adjust  the  duty  cycle  of  the  switch 
in  the  converter  [25-27].  This  section  describes  the  design  of  the 
compensation  networks. 

Assume  that  all  power  switches  and  diodes  are  ideal  elements, 
and  the  circuit  operates  in  continuous-conduction  mode  (CCM). 
Additionally,  an  equivalent  series  resistance  rc  associated  with  the 
capacitor  and  a  series  parasitic  resistance  rL  associated  with  the 
inductor  including  the  converter  equivalent  circuit  as  shown  in 
Fig.  2  are  added  to  increase  the  simulation  accuracy  [28-30].  As 
seen  from  Fig.  2,  the  minimum  values  of  the  inductances  and  capac¬ 
itances  in  the  converter  can  be  obtained  from  the  input  voltage, 
load  size,  switching  frequency,  and  output  voltage.  These  values 
are  used  as  a  reference  to  design  the  electrical  elements  of  power 
converters.  Since  the  converter  itself  is  a  part  of  the  control  loop,  the 
design  of  such  control  processes  requires  a  knowledge  of  the  small- 
signal  characteristics  of  the  converter.  Based  on  the  gain  and  phase 
margin  of  the  open-loop  control-to-output  transfer  function  of  the 
converter,  the  compensator  is  then  utilized  to  mitigate  steady-state 


Table  1 

Parameters  of  the  Ballard-Mark-V  PEMFC. 


Symbol 

Parameter 

Value 

Unit 

^an.in 

Upstream  nozzle  coefficient 
of  a  node  channel 

4.2  x  10“5 

kg  (s  atm)-1 

kan.out 

Downstream  nozzle 
coefficient  of  a  node  channel 

4.2  x  10-6 

kg  (s  atm)-1 

Pan,s 

Pressure  of  anode 

35 

psig 

T 

1  an, in 

Upstream  hydrogen 
temperature  of  anode 
channel 

296.5 

K 

cv,H2 

Specific  heat  of  hydrogen  at 
constant  volume 

10,124.71 

J(kgK)-1 

CP ,  2 

Specific  heat  of  hydrogen  at 
constant  pressure 

14,209 

J(kgK)-1 

Van 

Volume  of  the  anode  channel 

0.005 

m3 

k  conv.an 

Heat  convection  coefficient 
between  the  cell  body  and 
the  anode  channel 

2 

WK-1 

Rh2 

Gas  constant  of  hydrogen 

4071.1  x  10-5 

atm  m3  (kgK)-1 

Mh2 

Molecular  mass  of  hydrogen 

2  x  10-3 

kg  mol-1 

^ca,in 

Upstream  nozzle  coefficient 
of  a  cathode  channel 

1.8  x  10-3 

kg  (s  atm)-1 

kca.out 

Downstream  nozzle 
coefficient  of  a  cathode 
channel 

2.71  x  10“4 

kg  (s  atm)-1 

Pc  a,s 

Pressure  of  cathode 

35 

psig 

T 

1  ca,in 

Upstream  hydrogen 
temperature  of  cathode 
channel 

296.5 

K 

cv,02 

Specific  heat  of  oxygen  at 
constant  volume 

662.2 

J(kgK)-1 

cp,02 

Specific  heat  of  oxygen  at 
constant  pressure 

922 

J(kgK)-1 

Vca 

Volume  of  the  cathode 
channel 

0.01 

m3 

kconv.ca 

Heat  convection  coefficient 
between  the  cell  body  and 
the  cathode  channel 

10 

WK-1 

Ro2 

Gas  constant  of  oxygen 

262.68  x  10-5 

atm m3  (kgK)-1 

Mo2 

Molecular  mass  of  oxygen 

32  x 10-3 

kg  mol-1 

^conv.amb 

Heat  convection  coefficient 
between  the  cell  body  and 
the  ambient 

17 

WK-1 

ft1bodyCp,body 

Product  of  mass  and  the 
average  specific  heat  of  the 
cell  body 

3.5  x  104 

JK-1 

Patm 

Ambient  pressure 

1 

atm 

famb 

Ambient  temperature 

296.5 

K 

cp,h20 

Specific  hest  of  water  at 
constant  pressure 

2186 

J(kgK)-’ 

A  HRj 

Lower  heating  value 

1.196  x  108 

J 

X 

Membrane  resistivity 
parameter 

12.5 

Coefficients  of  activation 
over  voltage 

-0.948 

_ 

£ 3 

Coefficients  of  activation 
over  voltage 

-1.93  x  10-4 

_ 

$4 

Coefficients  of  activation 
over  voltage 

7.6  x  10-5 

_ 

Qi 

Double  layer  capacitance 

0.035  x  232 

F 

hnem 

Membrane  thickness 

1.78  x  10-2 

cm 

E  o 

Reference  potential  at  unity 
activity 

1.229 

V 

R 

University  gas  constant 

8.314 

J(molK)-1 

hlcell 

Number  of  cells  in  the 
fuel-cell  stack 

35 

— 

-^cell 

Fuel  cell  active  area 

232 

cm2 

F 

Faraday  constant 

96,485 

Cmol-1 

errors.  As  a  result,  the  stable  output  voltage  can  be  obtained  from 
this  process. 

3.1.  Procedure  and  criterion  of  compensator  design 

The  procedure  of  compensator  design  is  simply  described  as 
follows: 
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Fig.  2.  Schematic  diagram  of  the  designed  buck  or  boost  converter. 


(1 )  simulating  the  frequency  response  of  the  open-loop  control-to- 
output  transfer  function  of  the  converter, 

(2)  setting  the  gain  crossover  frequency  for  the  converters  and 
using  this  frequency  to  feedback  the  gain  of  the  compensation 
circuit, 

(3)  selecting  the  compensator  type  as  shown  in  Fig.  2  and 
designing  the  electrical  elements  of  the  compensation  circuit, 
and 

(4)  simulating  small-signal  analysis  of  the  frequency  response  of 
the  closed-loop  control-to-output  transfer  function  of  the  con¬ 
verter. 

As  for  the  criterion  of  compensator  design,  it  includes  the  fol¬ 
lowing. 

(1 )  The  large  dc  gain  is  required  to  decrease  the  steady  state  error. 

(2)  The  small  gain  switching  frequency  is  required  to  suppress  the 
noise  produced  during  the  switching. 

(3)  The  gain  crossover  frequency  of  the  converter  is  selected  to  be 
between  1/6  and  1/10  of  the  gain  switching  frequency. 

(4)  The  phase  and  gain  margins  are  suggested  to  be  45-60°  and 
6-20  dB,  respectively. 

3.2.  Frequency  response  analysis 

As  for  the  designed  buck  converters,  the  slope  of  the  frequency 

response  of  the  open-loop  control-to-output  transfer  function  of 


the  converter  at  the  first  corner  is  -40dBdec_1  because  the  two 
poles  generated  from  the  LC  filter  fall  close  to  each  other.  Moreover, 
due  to  the  equivalent  series  resistance  associated  with  a  capaci¬ 
tor,  a  zero  occurs  there  at  high  frequency,  resulting  in  the  slope 
of  -20dBdec_1.  As  for  the  designed  boost  converter,  it  not  only 
has  the  characteristics  of  the  frequency  response  of  the  designed 
buck  converter  but  also  has  a  right-half-plane  zero.  Based  on  the 
Bode  diagram  (i.e.,  gain  and  phase  plots)  for  the  open-loop  control- 
to-output  transfer  function,  one  can  select  the  compensator  type 
as  shown  in  Fig.  2,  corner  frequency,  and  gain  curve.  The  values 
of  electrical  elements  used  in  the  compensation  circuit  are  listed 
in  Table  2  and  the  close-loop  simulation  circuit  with  the  com¬ 
pensator  for  the  designed  buck  and  boost  converters  is  shown  in 
Fig.  2. 

4.  SIMULINK  block  diagrams 

The  mathematical  expressions  for  the  fuel-cell  system  and 
power  converters  were  simulated  in  the  Matlab/SIMULINK™  envi¬ 
ronment.  A  block  diagram  for  building  a  dynamic  model  of  PEMFC 
and  power  converters  is  given  in  Fig.  3.  In  this  block  diagram,  the 
input  quantities  are  anode  and  cathode  pressures  (Pan  and  Pca)  and 
operating  temperature  of  the  fuel  cell  (Top).  Additionally,  mass  dif¬ 
fusion  equations  are  used  to  calculate  the  effective  partial  pressures 
of  H2  and  02.  Then,  the  Nernst  equation  is  employed  to  determine 
the  internal  potential  of  the  fuel  cell.  The  ohmic  voltage  drop  and 
activation  voltage  drop  together  with  the  double-layer  charging 
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Table  2 

Values  of  electrical  elements  in  the  circuits  of  the  compensators  and  converters. 


Parameters 

Converters  with  the  output  voltage  of 

48  V 

24  V 

12V 

5  V 

l^stack 

20-40 V 

24-40 V 

20-40 V 

20-40 V 

fs 

50  kHz 

50  kHz 

50  kHz 

50  kHz 

L 

10  |mH 

37.5  jjlH 

20  jjlH 

230  jjlH 

rL 

10  m^ 

10  m^ 

10  m^ 

10  m^ 

C 

1.6  mF 

70  |jlF 

600  |xF 

100[xF 

rc 

10m£2 

10mS2 

10  mQ 

10  mQ 

R\oad 

8.2-1.37  £2 

24-3.84  Q 

2.4-0.384  Q 

25-4  Q 

37.3  kQ 

14.5  k^ 

3.55  kQ 

37.8  kQ, 

R2 

316 

21.4  kQ 

15  k£2 

69.5  kQ 

R3 

1  kQ 

100  Q 

100^ 

W0Q 

Cl 

1.88  nF 

4.8  nF 

61.7  nF 

8  nF 

c2 

168  pF 

7  nF 

14.6  nF 

4.36  nF 

c3 

15.6nF 

- 

- 

- 

effect  are  applied  to  determine  the  output  voltage  of  a  single  fuel 
cell. 

In  Matlab/SIMULINK™  a  set  of  equations  can  be  grouped  and 
are  called  subsystems.  Detailed  of  these  subsystems  are  given  in 
Figs.  4-7.  The  PEMFC  subsystem  using  the  block  diagram  men¬ 
tioned  in  Fig.  3  is  shown  in  Fig.  4.  The  subsystem  of  reactant 
gases  diffusion  in  the  electrodes  given  in  Fig.  5  employs  (28)-(37) 
and  calculates  hydrogen/oxygen  partial  pressures.  The  subsystem 
of  electrochemistry  for  Ballard-Mark-V  PEMFC  employs  (38)-(48) 
using  the  parameters  in  Table  2  as  shown  in  Fig.  6.  It  takes  hydro¬ 
gen/oxygen  partial  pressures,  cell  current,  and  temperatures  as 
inputs  and  calculates  the  cell  voltage.  A  gain  at  the  terminal  of  the 
subsystem  multiplies  this  output  by  the  number  of  cells  inside  one 
stack.  The  subsystem  of  mass  and  energy  transfer  at  the  anode, 
cathode  and  cell  body  for  Ballard-Mark-V  PEMFC  using  the  control 
volume  approach  employs  (12),  (19)  and  (27)  with  the  parameters 
listed  in  Table  1  as  shown  in  Fig.  7. 

As  for  the  power  converter  subsystem,  it  includes  the  circuits 
of  power  converters  and  compensators  with  voltage  feedback  con¬ 
trol.  Based  on  the  electrical  loads  of  the  AUV,  four  converters  are 
designed  to  obtain  5-,  12-,  24-  and  48-Vdc  output,  and  compen¬ 
sators  are  designed  to  regulate  these  dc  voltages  in  the  presence  of 
transients.  In  Matlab/SIMULINK™  a  model  of  the  compensator  cir¬ 
cuit  with  voltage  feedback  control  is  built  based  on  the  references 
[31,32].  To  limit  the  paper  length,  Fig.  8  only  shows  the  model  of 
the  output  voltage  of  24V  dc-dc  converter.  As  seen  from  Fig.  8,  the 
repeating  sequence  block  is  used  to  generate  the  sawtooth  wave¬ 
form,  and  the  relay  in  discontinuities  blockset  is  used  to  simulate 
the  comparator.  If  the  difference  between  the  compensating  signal 
and  the  sawtooth  signal  is  positive,  the  switch  is  open  and  the  out¬ 
put  value  is  one;  if  the  difference  between  the  compensating  signal 


and  the  sawtooth  signal  is  negative,  the  switch  is  closed  and  the 
output  value  is  zero. 

5.  Simulations  results 

5.2.  Effects  of  operating  conditions  on  cell  performance  of  PEMFC 

Fig.  9  shows  the  effects  of  using  air  or  pure  oxygen  as  the  oxi¬ 
dant  on  the  output  voltage  and  output  power  for  Ballard-Mark-V 
PEMFC.  As  seen  from  Fig.  9,  the  output  voltage  and  output  power 
using  pure  oxygen  as  the  oxidant  are  higher  than  those  obtained  for 
the  air.  Due  to  the  low  operating  temperature  of  the  Ballard-Mark- 
V  PEMFC,  an  increase  in  the  operating  temperature  may  result  in  a 
better  cell  performance  because  of  a  higher  efficiency  of  the  electro¬ 
chemical  reaction.  Fig.  10  shows  the  output  voltages  of  the  studied 
fuel  cell  under  different  operating  conditions  of  temperature  and 
pressure.  As  seen  from  Fig.  1 0,  a  linear  increase  in  the  output  voltage 
with  increasing  the  operating  temperature  and  pressure  has  been 
observed  in  a  high  current  density  range.  In  general,  the  operat¬ 
ing  temperature  of  the  studied  fuel  cell  is  limited  to  approximately 
80  °C,  and  both  the  anode  and  the  cathode  sides  operating  pres¬ 
sures  are  from  1  to  3  atm.  As  a  result,  due  to  the  limitations  of 
the  underwater  environment,  the  PEMFC  uses  pure  hydrogen  and 
pure  oxygen  as  fuel  in  this  study.  Assuming  the  pressures  of  oxy¬ 
gen  and  hydrogen  in  the  inlet  gas  are  operated  at  3  atm  with  the 
temperature  at  approximately  70  °C. 

5.2.  Simulation  of  PEMFC  subsystem 

Fig.  1 1  shows  the  total  output  voltage  and  output  power  for 
Ballard-Mark-V 35-cell  5-kW  PEMFC  based  on  the  parameters  men¬ 
tioned  in  Table  2.  As  seen  from  Fig.  1 1 ,  the  sharp  voltage  drop  at  low 
currents  is  caused  by  an  activation  overvoltage,  and  the  linear  drop 
at  intermediate  currents  is  due  to  an  ohmic  overvoltage.  Moreover, 
the  maximum  output  power  occurring  at  the  load  current  of  239  A 
is  about  3943  W,  and  the  corresponding  output  voltage  is  about 
16.5  V. 

Since  the  loads  of  the  designed  AUV  as  listed  in  Table  3 
are  below  21<W,  the  maximum  output  current  of  the  studied 
fuel  cell  is  computed  to  be  about  78  A,  and  its  corresponding 
output  is  ranged  from  27  V  to  40  V.  Fig.  12  shows  the  out¬ 
put  voltage  and  output  current  of  Ballard-Mark-V  35-cell  5  kW 
PEMFC  at  a  power  output  between  0.3  kW  and  2kW.  As  seen 
from  Fig.  12,  the  PEMFC  output  voltage  is  not  stable  due  to 
the  polarization  losses  (i.e.,  activation  and  ohmic  overvoltages). 
Thus,  PWM  dc-dc  converters  (i.e.,  three  buck  converters  and  one 
boost  converter)  are  connected  to  the  PEMFC  output  for  con¬ 
verting  unregulated  dc  voltage  into  four  desired  voltage  levels. 
Note  that  PWM  dc-dc  converters  are  designed  to  operate  under 
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Pea 

Pan 


Fig.  3.  Diagram  of  building  the  models  of  a  PEMFC  and  dc-dc  power  converters. 
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Fig.  4.  SIMULINK  block  diagram  of  the  Ballard-Mark-V  PEMFC  subsystem. 


overloads  and  light  load.  The  light-load,  full-load,  and  overload  con¬ 
ditions  are  defined  to  be  20%,  100%  and  120%  of  the  maximum 
rated  power  of  the  designed  AUV  as  listed  in  Table  4,  respec¬ 
tively. 

For  the  purpose  of  simulation,  no  loss  of  energy  is  assumed  for 
converting  an  input  energy  of  Ballard-Mark-V  35-cell  5-kW  PEMFC 
into  an  output  power.  The  output  power  and  output  voltage  of  the 
studied  fuel  cell  during  a  load  change  from  full  load  to  light  load  for 
all  four  buses  of  voltages  are  shown  in  Figs.  13  and  14.  As  seen  from 
Fig.  13,  the  total  output  power  of  the  studied  fuel  cell  varies  from 
1 841  W  to  370  W  as  a  load  change  occurs  at  4  s  during  a  simulation 


time  of  8  s.  As  seen  from  Fig.  14,  the  output  voltage  and  output 
current  of  the  studied  fuel  cell  are  29.3  V  and  63.1  A  at  full  load 
condition,  respectively.  However,  the  output  voltage  and  output 
current  of  the  studied  fuel  cell  are  32.5  V  and  11. 4  A  at  light  load 
condition,  respectively. 

5.3.  Simulation  ofPWM  dc-dc  power  converter  subsystem 

As  for  the  designed  AUV,  four  PWM  dc-dc  converters  with  the 
compensation  network  using  voltage-model  control  are  connected 
to  the  output  of  the  studied  fuel  cell  for  converting  unregulated 


Fig.  6.  SIMULINK  block  diagram  of  the  electrochemistry  subsystem  of  Ballard-Mark-V  PEMFC. 
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Fig.  7.  SIMULINK  block  diagram  of  mass  and  energy  transfer  for  the  Ballard-Mark-V  PEMFC. 


voltage  into  four  desired  voltage  levels  during  a  load  change.  Gen¬ 
erally,  the  converter  output  voltage  will  dip  at  the  sudden  change  of 
load  from  light  to  heavy  or  jump  at  the  sudden  change  of  load  from 
heavy  to  light.  The  magnitude  and  duration  of  the  transient  voltage 
occurring  in  the  output  voltages  during  a  load  change  are  designed 
to  be  within  ±5%  and  3  ms  in  this  paper.  However,  the  range  of 
acceptable  variation  in  the  output  voltages  usually  depends  upon 
the  users’  applications.  The  output  voltages  and  output  currents 


for  5-,  12-,  24-  and  48-V  dc-dc  converters  during  a  load  change  are 
shown  in  Figs.  15-18,  respectively.  Since  the  output  of  the  studied 
fuel  cell  ranges  from  27  V  to  40  V,  only  the  48-V  dc-dc  converter  will 
need  to  use  a  boost-converter.  The  rest  of  converters  will  use  buck 
converters.  As  seen  from  Figs.  15-18,  the  output  voltages  of  four 
dc-dc  converters  all  have  a  transient  spike  during  a  load  change  at 
4  s.  Nevertheless,  their  output  voltages  are  all  stable  right  after  the 
load  change. 
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Fig.  9.  Output  voltage  and  output  power  of  the  Ballard-Mark-V  PEMFC  for  using  air 

or  pure  oxygen  as  the  oxidant  Fig.  H-  Output  power  of  the  Ballard-Mark-V  35-cell  5-kW  PEMFC  and  its  corre¬ 

sponding  output  voltage. 


5.4.  Simulation  of  the  supply  of  the  reactant  gases 

In  this  subsection,  a  parametric  analysis  to  simulate  transient 
dynamic  effects  on  the  supply  of  the  reactant  gases  within  the  fuel 
cell  system  is  carried  by  using  the  parameters  listed  in  Table  1. 
Based  on  Table  4,  the  output  power  of  the  fuel  cell  system  needs  to 


(a)  under  different  operating  temperatures  at  PH,  /  =  3  atm  /  3  atm. 


Table  3 

Voltage  levels  and  loads  for  the  designed  AUV. 


Rated  voltage 

Name  of  devices 

Maximum  rated  power 

DC  5V  Bus 

Embedded  control  main  board 

5  W 

DC  12V  Bus 

Inertial  measurement  unit 

0.78  W 

Global  position  system 

1  W 

Pressure  transducer 

0.24  W 

RF  modem 

9  W 

CTD  sensors 

0.252  W 

Headlamp 

1 00W  *  3 

DC  24V  Bus 

Ultra  short  baseline 

10W 

Acoustic  modem 

8  W 

Inertial  navigation  system 

6  W 

Doppler  velocity  log 

80  W 

Digit  video  camera 

5-10  W 

Forward-looking  sonar 

5-10  W 

DC  48V  Bus 

Thruster 

350  W  *  3 

Total  about  1841  W 

generate  about  1841  W  when  four  converters  are  operated  on  full 
load  condition,  and  it  only  needs  about  370  W  when  load  change 
from  full  load  to  light  load  is  simultaneously  occurred  on  four  con¬ 
verters.  The  time  length  for  the  simulation  is  900  s. 


c  50 
%  40 

? 30 

o> 

S  20 

> 


2.  10 


"5 

O  0 


- Output  voltage 

- Output  current 


_ r - V  ^ 


1 _ i _ i _ i _ »  1 

0  1  2  3  4  5  6 


Simulation  timo  (s) 


7 


'  200 

-  150 

-  100 
-  50 

:  o 

8 


Fig.  10.  Output  voltage  of  the  Ballard-Mark-V  PEMFC  under  different  operating  Fig.  12.  Current-voltage  curve  of  the  Ballard-Mark-V  35-cell  5-kW  PEMFC  at  the 
conditions.  power  output  between  0.3  and  20  kW. 


3820 


C.-H.  Lee,J.-T.  Yang / Journal  of  Power  Sources  196  (201 1 )  3810-3823 


Fig.  13.  Output  power  of  the  Ballard-Mark-V  35-cell  5  kW  PEMFC  when  the  load 
changes  from  full  load  to  light  load. 


Fig.  19  shows  the  temperature  histories  of  the  stack  body, 
anode  channel  and  cathode  channel.  The  temperatures  of  these 
three  control  volumes  exhibit  similar  trends.  During  the  first  300  s 
(i.e.,  at  light  load),  the  temperatures  slowly  increase.  The  rea¬ 
son  is  that  the  heat  generated  by  the  electrochemical  reaction  is 
more  than  the  heat  transferred  to  the  surrounding  atmosphere 
through  convection.  During  the  second  stage  (300-600  s)  (i.e.,  load 
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Fig.  14.  Current-voltage  curve  of  the  Ballard-Mark-V  35-cell  5  kW  PEMFC  when  the 
load  changes  from  full  load  to  light  load. 


Table  4 

Classification  of  light  load,  full  load  and  overload  of  the  designed  AUV  at  each  dc 
bus. 


Rated  voltage 

Light  load  (20%) 

Full  load  (100%) 

Overload  (120%) 

DC  5V 

1  W 

5  W 

6.25  W 

DC  12V 

60  W 

312  W 

375  W 

DC  24V 

24  W 

124  W 

150  W 

DC  48V 

280  W 

1400  W 

1680  W 
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Fig.  15.  Output  voltage  and  output  current  of  5  Vdc  converter  during  a  load  change  from  5  £2  to  25  £2. 
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Fig.  16.  Output  voltage  and  output  current  of  12  Vdc  converter  during  a  load  change  from  0.46  £2  to  2.4  £2. 
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Fig.  17.  Output  voltage  and  output  current  of  24  Vdc  converter  during  a  load  change  from  0.46  £2  to  24  £1 
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Fig.  18.  Output  voltage  and  output  current  of  48  Vdc  converter  during  a  load  change  from  1.46  £2  to  8.2  £2. 


change  from  light  to  full  load),  the  output  power  of  the  stack 
increases  due  to  the  load  current  increase,  and  correspondingly 
the  heat  generated  by  the  electrochemical  reaction  increases  and 
overcomes  the  heat  loss  due  to  convection.  Thus,  the  tempera¬ 
tures  of  stack/anode/cathode  all  increase.  During  the  last  stage 
(600-900  s)  (i.e.,  load  change  from  full  to  light  load),  the  temper¬ 
atures  drop,  as  the  output  power  and  hence  the  heat  generated 
is  reduced.  It  is  worth  mentioning  that  the  temperatures  of  the 
three  control  volumes  differ  significantly.  This  clearly  demonstrates 
the  advantage  of  dividing  the  fuel  cell  system  into  three  control 
volumes. 

Fig.  20  shows  the  pressure  variation  within  the  two  chan¬ 
nels.  When  the  load  current  is  high  (i.e.,  load  change  from  light 
to  full  load),  the  fuel  consumed  by  the  electrochemical  reaction 
increases.  Correspondingly,  the  gas  pressures  go  down.  On  the  con¬ 
trary,  when  the  load  current  is  low  (i.e.,  load  change  from  full 
to  light  load),  the  pressures  within  the  channels  become  rela¬ 
tively  high.  Fig.  21  shows  the  change  of  mass  of  gases  (hydrogen, 
oxygen  and  nitrogen)  within  the  two  channels  due  to  the  load 
current  change.  During  the  first  stage  and  last  stage,  the  load  cur¬ 


rent  is  lower,  which  leads  to  lower  fuel/oxygen  consumption  for 
the  electrochemical  reaction  and  higher  mass  of  gases  within  the 
channels.  The  opposite  can  be  observed  for  the  second  stage.  It 
is  shown  that  hydrogen  and  oxygen  decrease  due  to  load  current 
increase.  This  is  because  when  the  oxygen  consumption  increases, 
the  pressure  within  the  cathode  channel  goes  down  and  more 
oxygen  will  enter  from  the  inlet.  The  inlet  mass  flow  rates  of 
hydrogen  and  oxygen  are  shown  in  Fig.  22.  Correspondingly,  their 
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Fig.  19.  Temperature  change  at  the  anode  channel,  cathode  channel  and  cell  body. 
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Fig.  20.  Pressure  change  at  the  anode  and  cathode  channel. 
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Fig.  21.  Hydrogen  and  oxygen  change  within  the  anode  and  cathode  channels. 


(b) 

Fig.  23.  Mass  consumption  of  hydrogen  and  oxygen. 


6.  Conclusions 

This  paper  studies  the  electrical  behavior  of  Ballard-Mark-V  35- 
cell  5  kW  PEMFC  with  its  output  voltage  range  from  27-  to  40-V 
based  PWM  dc-dc  power  conversion  under  load  changes  for  appli¬ 
cations  in  AUVs.  The  models  of  the  studied  fuel  cell  and  PWM  dc-dc 
converters  have  been  implemented  in  Matlab/SIMULINK™  envi¬ 
ronment.  The  PEMFC  model  includes  the  double-layer  charging 
effect,  gas  diffusion  in  the  electrodes,  and  the  thermodynamic  char¬ 
acteristic  while  the  operation  of  PWM  dc-dc  converter  assumes  to 
be  in  CCM  with  a  compensator  in  voltage-mode  control.  A  series 
of  simulations  show  that  the  output  voltages-based  dc-dc  power 
conversion  with  the  designed  compensation  network  is  stable  as 
the  load  changes  from  full  load  to  light  load.  However,  we  regret 
that  the  results  of  the  simulations  have  not  yet  been  compared  with 
actual  measurements  since  the  experimental  work  is  in  progress. 


mass  consumptions  are  shown  in  Fig.  23.  The  simulation  results 
presented  here  clearly  illustrates  the  complicated  dynamic  inter¬ 
actions  between  temperature  and  gas  flow,  and  demonstrates  the 
necessity  of  developing  a  transient  model  using  the  control  volume 
approach. 


(b) 

Fig.  22.  Inlet  mass  flow  rate  of  hydrogen  and  oxygen. 
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